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Abstract 

An analytic resolution of the covariant oscillator quark model for a three-body system is pre- 
I sented. Our harmonic potential is a general quadratic potential which can simulate both a A-shape 

o ■ 

^Sj ' configuration or a simplified Y-configuration where the junction is located at the center of mass. 
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The mass formulas obtained are used to compute glueball and baryon spectra. We show that 
the agreement with lattice and experimental data is correct if the Casimir scaling hypothesis is 
assumed. It is also argued that our model is compatible with pomeron and odderon approaches. 



(N 

(3 ; PACS numbers: 12.39.Pn, 12.39.Ki, 12.39.Mk, 14.20.-c 

(N 



in 
o 



*FNRS Research Fellow; E-mail: "fabien .buisseret @umh. a c .bel 
^^FNRS Research Associate; E-mail: claude . semay @ umh . ac . be| 

Typeset by REVT^ 1 



I. INTRODUCTION 



The covariant oscillator quark model (COQM) is a phenomenological model of hadrons 
It is not based on a Hamiltonian like usual quark models, but on an operator giving the 
square mass of the considered system. This model is not a QCD inspired one like a spinless 
Salpeter Hamiltonian with a linear confining potential. But, through a covariant formalism, 
it is able to deal in a simple way with retardation effects, which are rather complicated 
in other approaches . Moreover, the COQM gives the correct Regge slopes of the 

,„eso„,c a„d baryo„,c Regge trajectodes M. An attempt to reproduce the ,„ass spectrum 
of glueballs formed of two gluons with this model has also been made in Ref. [2|- The 
key ingredient of the COQM is the particular form of the confining potential. In a two- 
body system, it is a harmonic potential depending on the separation of the particles, the 
spring constant of the potential determining the value of the Regge slope. For the baryons, 
a quadratic equivalent of the so-called A-potential has been used in Ref. joj. However, 
recent developments in lattice QCD rather support the picture of an Y-junction inside the 
baryons j^, and thus the potential of the three-body COQM should be modified. Moreover, 
accurate studies have since been performed in lattice QCD about the spectrum of glueballs 
0, [l^. The purpose of this paper is thus to reconsider the COQM taking into account 
these new results. To do this, we will use a more physically relevant three-body potential, 
simulating an Y-junction, and we will show that the mass spectra and the wave functions 
can be analytically found. This is somewhat unusual in a three-body problem. Solving the 
equations of the COQM will give us the possibility to compare its predictions with the well 
known baryon spectra, and new lattice QCD results concerning the glueballs. Moreover, it 
will give us a first estimation of the retardation effects in a three hadrons system, which we 
hope will open the way for a three body generalization of the work in Ref. ^ \ . 

Our paper is organized as follows. In Sec. El we introduce the formalism of the COQM 
through the simple case of a two-body system, as it is done in Ref. Then, in Sec. IIIH 
we discuss the introduction of a general potential for the three-body COQM, which can be 
seen as a mixing of a harmonic equivalent of the Y-junction and the harmonic potential of 
Ref. Q]. After that, we show that it is possible to find analytical solutions for the three-body 
COQM with our new potential. In Sec. IIVI we discuss the value of the parameters of our 
model and we compare its results with the experimental and lattice QCD data. Finally, we 
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sum up our results in Sec. El 



II. TWO-BODY COQM 

The equation of the COQM for interacting hadrons is 



k=l 



Pk 

2mk 



0, 



(1) 



where the potential U, assumed to depend only on relative coordinates, describes the confin- 
ing interaction between particles whose four- vector coordinates are x^. We can perform 
a change of coordinates and express the Xi in terms of the center of mass position, R, and 
— 1 relative coordinates, r^. Then, Eq. (P) reads 



7ri 



2mj 



a=l ^^^^ 



0, 



(2) 



where P and tt^ are the momenta associated to R and respectively. The /Iq are the 
reduced masses and rrit = nii + ■ ■ ■ + m„. But, we know that = — M^, where M is the 
mass of the system. So, instead of a Hamiltonian, the COQM allows us to write an equation 
giving the square mass of the system from Eq. ^ 
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(3) 



In this section, we will only consider the quark-antiquark system. The confining potential 
for mesons is ^ 

U = ^K^r^ - W, (4) 

with r = Xi — X2 = (a, r, 6, ip). The spring constant K2 is a parameter of the model, as well 
as the constant W , which can take into account in a very simple way other contributions 



than the confinement: One gluon exchange, spin interactions,. . . 
of Eq. (jni) can be written in this CcLSG cLS 
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The quantized version 



2mt 



2/i 2 



W 



(5) 



with /i = 17111712 /iTLt the reduced mass. Applying the well known theory of the harmonic 
oscillator, one finds 

M"^ = 2mt^ICJJl{2n + ^+l)-2mtW, (6) 



^ = Xo(o-)0n/r)Fr(^,^), (7) 



where the are the spherical harmonics. With (3 = v^/i^2, 

1/4 



is an eigenfunction of the one dimensional harmonic oscillator, and 



is a radial eigenfunction of the three dimensional harmonic oscillator [12] • Hn and L" are 
the Hermite and Laguerre polynomials respectively. It is worth mentioning that the usual 
factor (2n + £ + 3/2) of the harmonic oscillator is here replaced by (2n + £ + l) in Eq. To 
understand this, we have to consider the contribution of the relative time a. When Eq. (0) 
is solved, two harmonic oscillators appear: one for the spatial part and one for the relative 
time part. The second has the opposite sign of the first and will contribute to decrease the 
square mass. But a nonphysical degree of freedom is now present: which eigenstate of the 
relative time oscillator do we have to choose? The prescription, which can be written in a 
covariant way, is to consider only the fundamental state as a physical one That is 

why we find Xo in relation ((Tj) and —mt\jK2/fi for the contribution of the relative time. 
This contribution is in fact due to the retardation effects as they appear in the COQM. The 
most probable value for a is 0, which is in agreement with the hypothesis of Ref. When 
the two particles have the same mass m, formula (jHl) reduces to 



= ^J?>2mK2 (2n + £ + 1) - 2mW. (10) 

It has been shown in Ref. ^ that the COQM is able to predict the meson Regge tra- 
jectories in agreement with the experimental data with K2 = 0.107 GeV^. We see from 
Eq. ® that the COQM should not be used for heavy mesons, because only light mesons 
exhibit Regge trajectories. Finally, we can remark that the COQM is not relevant for mass- 
less particles. This drawback is inherent to this approach, but it is not troublesome if the 
constituent quark and gluon masses are used instead of the current ones. 
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III. THREE-BODY COQM 



A. The confining potential 

As in the previous section, we will only consider quark systems, that is to say baryons. 
We summarize here some considerations of Ref. The potential which is used in this 
paper has the form 



It can be seen as an harmonic equivalent of the usual A-potential. With the same quark 
masses as in the mesonic case, a good agreement between theoretical and experimental 
baryonic Regge slopes can be obtained if = The proposed justification of 

such a 1/4 factor is the following heuristic argument: in a usual meson model with linear 
confinement, the potential term of the Hamiltonian reads a2r, where 02 is the energy density 
of the flux tube between the quark and the antiquark and r is their spatial separation. In the 
COQM, the potential appearing in the square mass operator is K2r^- So, one can consider 
that K2 oc a\. Then, if we assume that the energy density 03 of the flux tube between two 
quarks i and j in a baryon is proportional to the color Casimir operator AjAj/4, we should 
have 03 = 02/2 and thus K^^ = -ft'2/4 for a A-shape potential. 

However, recent developments in lattice QCD tend to confirm the Y-junction as the 
more realistic configuration for the flux tube in the baryons 8]. In this picture, flux tubes 
start from each quark and meet at the Toricelli point of the triangle formed by the three 
particles. This point, denoted xt, is such that it minimizes the sum of the flux tube lengths, 
and its position is a complicated function of the quark coordinates Xj. Moreover, the energy 
density of the tubes appears to be equal for mesons and hadrons: 03 = 02. As we want to 
include these recent results in the COQM, we have to change the expression of the harmonic 
potential. In particular, we have to keep = K2 and to take a quadratic equivalent of the 
Y-junction. 

In Ref. the complex Y-junction potential 

3 

= 03^ |fi - xtI (12) 

i=l 




(11) 
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is approximated by the more easily computable expression 



1=1 



R 



(13) 



where R is the position of the center of mass. If a = 1, Eq. ()13|1 is a simplified Y-junction, 
where the Toricelli point is replaced by the center of mass. If a = 0, this interaction reduces 
to a A-potential. Let us note that the presence of the factor 1/2 in the A-part of the 
potential is purely geometrical and simply arises because in a triangle ABC with a Toricelli 
point T, \AT\ + \BT\ + \CT\ > {\AB\ + \BC\ + \CA\)/2. Results of Ref. jisj, obtained in 
the framework of a potential model, show that a = 1 gives better results than a = 0, and 
that the Y-junction is approximated at best when a is close to 1/2. 

In order to simulate at best the genuine Y-junction, keeping the calculations feasible, 
we assume, in agreement with Eq. (|TH|l . the following expression for the potential in the 
baryonic COQM 



2 



a ^ ( Xi 

i=l 



Rf 



a 



Kj 



(14) 



with K3 = K2. The origin of the factor 1/4 = (1/2)^ is now seen as geometrical only. We 
define 

= aKs, p = Ksil - a)/4, (15) 

and the potential (fT^ becomes 



U ^^2_^(x,- R) 

i=l ^ i<j 

which is the expression we use in the following. 



Xi 



Xi 



(16) 



B. Mass formula and wave function 



We have to solve Eq. (Q), which in this case reads 

3 



i=l 



with U given by Eq. (fT^ . First of all, we will replace the quark coordinates Xi = {xi, X2, x^} 

by x'l^ = {R, ^, T]}, with the center of mass defined as 

^ _ miXi + 1712X2 + m^x^ 
mt 
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rrit = mi + m2 + and 77} are two relative coordinates. The change of coordinates is 
made via a matrix Q, thanks to the relation Xi = Qn^x'f^. Let us note that the invariance of 
the Poisson brackets demands that Pi = iQ^^)JkP'k, with p[ = {P, 77^,77^}. We define 

^ 1 A 
1 C D 



Q 



(19) 



and find that the elements of Q can be constrained by the following equations 

A = -^C-^E, (20) 

nil mi 

B = -'^D-^F, (21) 
mi mi 

C = :^irni + mtED) , (22) 

D = FS, (23) 

mim2 [D{mi + m2) + Fma] , , 

-C/ = 2 T) [■^'*) 

mt [(Z}m2 + Fm^) + D'^mim2 + F'^mim^\ 

where 5 is a solution of 

mi - m2)m2{mtp + m-ip)5'^ + (m2 - ms) [m2msfi - mimt{fi + 2p)] 6 

+ {ms - mi)m3{mtp + m2p) = 0. (25) 

Constraints and (PT|) are consequences of the definition (fTHj) . Equation ensures that 
det Q = 1 in order to simplify the calculations of the p-. These three relations are sufficient to 
cancel the terms containing the cross products P-tt^ and P-tt^ in the kinetic part of relation 
p7|) . The last cross product tt^ • tt^ vanishes thanks to Eq. ()24|1 . Finally, formulas ()2H|1 and 
suppress the terms containing the cross product ^ • 77 in the potential (fT?)|) . 
The last parameter to fix is F, which has to be nonzero. We define 



[((5m2 + ms)^ + mi((5^m2 + ma)] /mi, (26) 

(27) 



mim2m^ 



mt 



and we choose 



F^ = 0/r. (28) 

We can then rewrite Eq. (|T7jl as 

P2 1 1 



+ TTT + 7^ + o^?^' + o^.^' - = 0, (29) 



2mt 20 20 2 2 



where 

= fi{A^ + C^ + E^)+p [{A - Cf + {A- Ef + {C- Ef] , (30) 



% = fi{B^ + D^ + F^)+p [{B - Df + {B- Ff + {D - Ff] . (31) 

Equation ^I^i has the nice property that the variables are all separated. Using the 
same arguments as those discussed in Sec. |Tll the square mass spectrum can be analytically 
computed. It reads 

= 2mt [{2n^ + + l)uj^ + {2n^ + ir^ + l)c^r, - W] , (32) 

with 



uj^ = y^^V^' f^i = ^^i^ ^r, = \/^r,/4>, Pn = 4>(^7j- (33) 
The internal wave function is given by 

where we used the definitions (jH} and (jHI). Only the ground state of the temporal part of 
the wave function must be considered, following the prescription of Ref. [^. It is worth 
mentioning that the most probable values for the relative times are C,^ = tj^ = 0, like in 
the mesons. It can be observed from Eq. ()32|) that the retardation effects bring a negative 
contribution to the squared mass, given by —mt{uj^ + cj^). 

As the COQM is expected only to be valid for light particles, the four possible baryonic 
systems are: nnn, snn, ssn or sss {n stands for u or d). Thus, we can always consider 
simplified cases where at least two masses are equal. When m2 = rris = m, a solution of 
Eq. ()25|1 is 5 = — 1, and one can find that 



, (mi + 2m)2p + (mf + 2m2)/i 

= \ , (35) 

" mimrrit 



_ /i + 3p 
'1 \l ■ 



(36) 



When the three masses are equal, we have simply 



= ^r, = \l = ^, (37) 
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and the square mass formula reduces to 

= 3v/(3 + a)K3m [2{n^ + n^) + {i^ + Q + 2] - 6mW. (38) 
It can be checked that for the particular case a = 0, our solution agrees with the one of 

n 

Ref. p. For three equal masses, a variation of a can be simply simulated by a variation of 
K3; this is not the case when particles have different masses. 



IV. COMPARISON WITH EXPERIMENTAL AND LATTICE QCD DATA 

Since the COQM includes neither the spin (S) nor the isospin (J) of the constituent 
particles, the data which can be reproduced here are the spin and isospin average masses, 
denoted M^v These are given by the following relation 

E„(2/+l)(2J+l)Af,,. 

E„(2/+1)(2JH-1) ' '-^^> 

with J = L + S , and where M/ j are the different masses of the states with the same orbital 
angular momentum i and the same quark content. We also use the formula (jH^ to compute 
a mass with the three-body COQM, but in this case / and J are replaced by the values of 
and irj corresponding to a given £ = £^ + 



A. Mesons and baryons 

It is a well known fact that the Regge slope of light mesons, such as nn states, is roughly 
equal to the Regge slope of the corresponding baryons nnn jlSl]. It is readily seen that the 
slope of the mass formula for meson (fTUI) and the slope of the mass formula for baryon ()38p 
are equal if 

a = ^ ^ 0.56, (40) 

a value close to the optimal one of 0.5 found in Ref. y^]. In the following, this value of a 
will be always used. 

Like in all potential models, the strength of the confining potential and the constituent 
quark masses must be fixed on data. We take the spring constants K2 = K3, as it is argued in 
Sec. mil and we fix the value of this parameter at 0.107 GeV^ as in Ref. [5|. The constituent 
quark masses are chosen to be equal to m„ = 0.313 GeV and nig = 0.375 GeV, in order 
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to obtain good baryon spectra (these values are different from those used in Ref. 0). The 
results are plotted in Figs. Hand|21 We can observe a good agreement between the COQM 
and the average square masses. What is important is that the correct slope is obtained, 
since the absolute values of the masses depend on ad hoc values of the parameter W. W is 
positive in every case, excepted for the sss baryons. This could be caused by strong spin 
interactions: masses of the sss states considered are not average ones. Our data concerning 
the baryons are taken from Ref. jlG^ . 

The difference between the constituent masses m„ and can seem too small, but such 
a mass difference is obtained in potential models in which a constituent state dependent 



mass is defined asTx/p'^ + z/^) where u is the current mass. The value of rru—mn is around 
40 MeV in Refs. |12l,llSl| for ground states of mesons and baryons. In Ref. an analytical 
approximate formula is given for the constituent quark mass of the baryon ground states : 
a difference of about 30 MeV is found for current masses 1/^ = and Ug = 130 MeV 3|. 
In our model, the constituent masses are not state-dependent but, in these references, small 
values for the difference rris — rrin are also found, and absolute masses m„ and are in 
agreement with our values. 

Formula (jH)) implies that the square meson masses have a linear dependence on the 
orbital angular momentum £ 

M^ = pi + Po. (41) 

This relation is also well verified experimentally. Using the data and the procedure of 
Ref. |20| , average square meson masses can be computed, and the relation (jlT|l can be fitted 
on these data to obtain an average experimental Regge slope f3exp.- With our parameters 
optimized for baryon spectra, a theoretical Regge slope /3coqm can be computed for mesons. 
It can be seen on Table Uthat the two slopes differ only by around 6% for various qq systems. 
By increasing the value of m„ by about 30 MeV, it is possible to improve the theoretical 
meson spectra. The price to pay is a slight deterioration of the baryon spectra. 



B. Glueballs 



We will assume here that the energy density of the flux tube starting from a particle 
i is proportional to the Casimir operator ^i,tI^''T^^\ with T^^^ being a 5'[/(3) generator 
in the corresponding representation. Several approaches tend to confirm this hypothesis 
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2lL l22j |. Then, if 02 is the energy density in a meson, we should have a2g = (9/4)a2 in a 
glueball formed of two gluons, and 03^ = a2g where a^g is the energy density in a three-gluon 
gluebaU. In the COQM, the corresponding spring constants K^g and K2g will thus be given 
by Ksg = K2g = (9/4)2/^2, where K2 is the spring constant for a meson. In order to simulate 
at best the Y-junction we will also use a = 5/9 as in the baryonic sector. 

The constituent gluon mass nig is a parameter of the model which must be fixed as the 
constituent quark masses. Because of the scarcity of reliable experimental data, we will 
determine it by using lattice-QCD results about two-gluon glueball spectrum.. We consider 
that all the positive charge conjugate states given in Refs. 0, Q are two-gluon glueballs. 
This in agreement with the results of the potential model of Ref. ^|. These states are 
shown in Table ITTl with the average square masses computed using Eq. (jHEI)- By fitting the 



mass formu^ 
usual ones 



a (ITUll on these data, we obtain nig = 0.770 ± 0.340 GeV, a value close to the 
Formula (fTUj) shows that the two-gluon glueball mass depends on 



22 



3Q- 



the product K2gnig. If we assume a color-dependence on the square- root of the Casimir 



operator 



25l ]. it is necessary to use a gluon mass around 1.7 GeV, which seems too heavy. 



According to the glueball-pomeron theory, the J'^^ two-gluon glueballs, having a max- 
imum intrinsic spin 5* coupled to a minimum possible orbital angular momentum £, stand 
on Regge trajectories. In Ref. the following relation is obtained 



J = 0.25 M, 



gg 



This result is close to the one obtained in Ref. |2 
noticeably different 



J = 0.36M,2,+ Q_^^ 



1.08. 

In Ref. |2 

0.80 



(42) 

the relation obtained is 
(43) 



The value of the ordinate at origin in this model depends on the value chosen for the strong 
coupling constant. If we use our average square two-gluon glueball masses (see Table Hlj) 
and if we assume J = i + 2 (maximum J-coupling) , we find 



J = (0.27 ± 0.06) M^g + (0.59 ± 0.49). 



(44) 



One can see that our slope is in agreement with of the slope found in Ref. l26l , but our 



ordinate at origin is more compatible with the corresponding one in Ref. |22| | 
Relation ()44|1 is in good agreement with the relation found in the COQM of Ref. 



Finally, 



J = 0.26 M;^ + 0.67. 



(45) 
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If the interaction between gluons was spin-independent, we could expect that the hghtest 
three-gluon glueballs are the states with a vanishing total orbital angular momentum and 
J^^ = 0^+, 1 , 3 [2^. The lattice calculations yield very different results This 
can be due to very large spin-orbit effects, which have been observed in two-gluon glueballs 
[2^. In this situation, we can expect that our spinless COQM can just describe the three- 
gluon glueballs with a vanishing total orbital angular momentum. An average square mass 
is obtained in Table ITTl with the 1 , 3 states of Ref. The ^ state of this work is 
assumed to be a two-gluon glueball (see above). 

The equivalent of the pomeron for three-gluon glueballs is called the odderon. J three- 
gluon glueballs, having a maximum intrinsic spin coupled to a minimum possible orbital 
angular momentum, are also expected to stand on Regge trajectories. The following results 
are predicted in Ref. j3l for a QCD effective Hamiltonian and for a non relativistic potential 
model respectively 

Jeff. = 0.23 M|^^- 0.88, (46) 
Jnr = 0.18 M|^^ + 0.25. (47) 

et us note that a slope for ggg states which is half the slope for gg states has been suggested 



3- 



If we fit our square mass formula with the unique data of Table [H] for three-gluon 
glueballs, and if we assume J = i + 3 (maximum J-coupling) we obtain 

J = (0.27 ± 0.06) M^gg - (1.57 ± 1.54). (48) 

The same slope is predicted for two-gluon and three-gluon glueballs, as in the case of quark 
systems, because we assume a = 5/9 and K^^g = K2g, and also because we choose for the 
gluon mass the value of 0.770 GeV found for the two-gluon glueballs. Our result is then 
compatible with the one from the QCD effective Hamiltonian of Ref. 3|- 
It is worth mentioning that the Regge trajectories presented in Refs. 
concern glueballs with given J^^ . Our model can only predict spin average square masses. 
So, the absolute masses given by our model cannot, in principle, be compared directly with 
masses of states with given J^*^. The ordinates at origin of Regge trajectories calculated 
with our model could then have a significant error, but we can expect a correct computation 
of the slope. 
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V. SUMMARY OF THE RESULTS 



In this paper, we have solved the covariant oscillator quark model applied to three-body 
systems, with a general quadratic potential. This interaction is a linear combination between 
a A-potential and a simplified Y-j unction where the Toricelli point is identified with the 
center of mass; it is believed to be a good approximation of the true Y-junction potential. 
We found analytic expressions for the mass spectrum and for the total wave function. We 
then have shown that our results are in quite good agreement with experimental and lattice 
data concerning baryons and glueballs provided the Casimir scaling is true. However, a 
more detailed study of the glueballs with the COQM requires to take into account the spin 
interactions, which are particularly strong in these particles. 
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FIG. 1: Average experimental square masses of the nnn states (empty circles) and sss states (full 
circles), compared with the theoretical results (straight lines). The nnn states are computed from 
the N and A baryons, and the sss states are the Q baryons iWnnn = 0.312 GeV and Wsss = —0.232 
GeV). 
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FIG. 2: Same as Fig. but for the nns (empty squares) and ssn states (full squares). The nns 
states are computed from the A and S baryons, and the ssn are the S baryons {Wnns = 0.379 GeV 
and Wssn = 0.223 GeV). 
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TABLE I: Regge slope (3 for mesons (see formula (|1T|)). /3exp. is the value obtained by averaging 
experimental data (see Ref. j2^) and /?coQM is our theoretical value (see formula The values 
of parameters W for the various systems are also given. 

State /3exp. (GeV-2) /5coqm (GeV-^) W (GeV) 

nn 1.13 1.04 0.370 

ns 1.16 1.09 0.220 

ss 1.19 1.13 0.065 



TABLE IL Average square masses, Mf^, for glueball states. The third and fourth columns show 
the different states used to compute the quantities Mf^. Data are taken from Ref. Q], except the 
mass of the 0* state which is taken from Ref. 9|l. 



State (n + 1)L 




M (GeV) 




99 IS 


0++ 


1.710 ±0.130 


5.183 ±0.667 




2++ 


2.390 ±0.150 




IP 


0"+ 


2.560 ±0.155 


8.762 ± 1.089 




2-+ 


3.040 ±0.190 




ID 


0*++ 


2.670 ±0.310 


12.567 ± 1.702 




3++ 


3.670 ± 0.230 




999 IS 


1— 


3.830 ± 0.230 


16.720 ± 1.967 




3~- 


4.200 ± 0.245 
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